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Formalism based on equilibrium statistical thermodynamics is applied to communication networks 
of decision making individuals. It is shown that in statistical ensembles for choice models, properly 
defined disutility can play the same role as energy in statistical mechanics. We demonstrate addi- 
tivity and extensivity of disutility and build three types of equilibrium statistical ensembles: the 
canonical, the grand canonical and the super-canonical. Using Boltzmann-like probability measure 
one reproduce the logit choice model. We also propose using g-distributions for temperature evo- 
lution of moments of stochastic variables. The formalism is applied to three network topologies of 
different degrees of symmetry, for which in many cases analytic results are obtained and numerical 
simulations are performed for all of them. Possible applications of the model to airline networks 
and its usefulness for practical support of economic decisions is pointed out. 
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I. INTRODUCTION 

Models used in demand analysis of transportation and 
communication problems usually assume that demand 
represents the result of decisions of each individual in 
the population. These decisions consist of choices made 
among finite sets of possibilities. As an example of a 
sequence of choices, in the context of airline transporta- 
tion demand, consider decision to be taken by potential 
air passenger. If the passenger wills to fly from the origin 
port O to destination port D, he has to choose among of- 
fers of carriers, accounting for many factors relevant for 
decision, as e.g. scheduled departure and arrival times 
and how they relate to his needs, flight durations, flight 
fares, numbers of stops and changes on the route (OD), 
availability of seats in first class, probability of delays, 
declared and expected quality of service and many other 
aspects, too numerous to itemize. In order to quantify 
such process of decision making, one incorporates dis- 
crete choice models in hope of better understanding and 
predicting behaviour of such complex system as trans- 
portation network and thus obtaining hints for marketing 
and revenue management decisions. 

One of the key concepts used in decision theory, dat- 
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ing back to the XlX-th century economy by Jeremy Ben- 
tham and other representatives of so called utilitarians 
economy, is the concept of utility, which had been fur- 
ther evolving since that time and found its more rigorous 
utterance in the classical work of von Neumann and Mor- 
genstern }M and in neoclassical economic theory (cf. e.g. 
refs 0, 0)- Utility function ascribes real numbers to 
alternatives and so defines the preference order on the 
choice set and puts the choice making process on more 
rational footing. 

An approach to communication based on discrete 
choice models has long history and abundant literature 
(cf. e.g. refs Q, Utility function in neoclassical 

economic theory is defined axiomatically. In order to use 
it in specific applications of choice models, one has to 
postulate a lot from outside of the theory and to find 
an effective way to either derive it from any fundamen- 
tal or general theory or to parametrize it and estimate 
from data. Derivation of the utility function from first 
principles represents major theoretical problem because 
of lack of fundamental theory of human behaviour. Nor- 
mally, one has to rely on partially justified models using 
power series approximations for utility functions, often 
restricted to linear or quadratic terms (cf. ref. [6(), and 
finding empirically relevant variables and estimators for 
parameters, if any data are available at all. Apart from 
those problems, numerous conceptual difficulties arise 
when applying utility functions to quantify the behaviour 
of the choice makers and to predict them. In particular, 
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classical Bentham's approach of maximizing the overall 
utility, integrated over individuals, when taking market 
decisions, is often a subject of serious objections Q. 

Being aware of a long-lasting debate around the mean- 
ing and the role of utility functions in economy, social sci- 
ences and, more generally, in game theory, in this paper 
we try to incorporate it into economic games on some- 
what different way. First, extending from the concept of 
utility, we find an analogy between the utility function, 
as it is used in choice models, and the energy function for 
physical systems. After specifying the system, consisting 
of an airline network and a set of passengers, and its state 
space, we discuss assumptions normally taken in theories 
of systems in thermal equilibrium and build quasistatic 
probability measures on the state space for equilibrium 
ensembles. This allows us to derive complete classical 
thermodynamics for the communication network. We 
consider the canonical and the grand canonical statistical 
ensembles for systems with fixed and variable numbers of 
clients. 

As a possible extension of the model, we propose to 
consider also the network topology and the number of 
network vertices as a multidimensional random variable 
and formulate generalized canonical ensemble where the 
existence of any connection between sites may fluctuate 
and is represented by a binomial random variable. In 
the frameworks of the canonical and the grand canoni- 
cal formalisms, we find some analytical results for ther- 
modynamics of particular network topologies and per- 
form numerical simulations. We study thermodynamical 
potentials, response functions and correlations between 
extensive variables in the system. The case of random 
topology is discussed in general terms only. It needs more 
study of a couple of specific and self-contained problems 
and this, altogether with numerical results, is relegated 
to future work. 

We also develop another approach, based on so called 
q-distributions or escort distributions, which allows us to 
monitor fluctuations and correlations of interesting char- 
acteristics of the network, their temperature evolution 
and phase relations in the network. 

Building our formalism we keep in mind its specific 
application to networks of airline traffic and often use 
its terminology, although we believe that this frame- 
work and most of its features are general enough to be 
used for other network-related problems, as various kinds 
of telecommunication, communications in social groups, 
traffic routing, energy networks, etc. On the other hand, 
the model we build does not account for many features of 
real airline markets, either because of lack of knowledge 
about utility functions or simplifications of the model at 
this stage. 



II. ASSUMPTIONS AND FUNDAMENTALS OF 
THE MODEL 

A. Axiomatics of utility 

In choice models considered here we assume that de- 
cision makers are individuals and that these individuals 
are independent of each other. Such models are usually 
called disaggregate. Assumption on independence of de- 
cisions between individuals is rather poorly justified and 
often just not true in reality. But for many practical ap- 
plications this assumption is not restrictive and depends 
on further details of the model, in particular on the util- 
ity function. In many cases, also in this paper, one may 
redefine the notion of a decision maker and consider a 
group of individuals as the decision maker, without loss 
of generality, and only mutual independence of such de- 
fined decision makers is relevant. This assumption can 
be further criticized because the ways decisions are taken 
by groups, even for very restricted class of consumer de- 
cisions considered here, certainly depend on the size and 
internal communication structure of groups and are qual- 
itatively different than in case of individuals. Thus the 
utility function for aggregated decision makers may differ 
from that of individuals. In order to give full account of 
these diffences in general case, which we do not pretend 
to do in this paper, one has to consider many additional 
effects, as internal negotiations inside groups, their dy- 
namics, possible splits, compactness of groups etc. How- 
ever, at this stage of model development and for the case 
of limited number of variables relevant for utility, and 
for reasonably homogeneous groups, we feel it is fair to 
assume that utilities for all decision makers are indepen- 
dent and identically distributed stochastic variables. In 
other words, we assume high degree of decision coher- 
ence of groups, in the sense that decisions of aggregated 
individuals do not differ significantly from that of real 
individuals. 

Analysis of the choice requires knowledge of the com- 
plete set of options disponsible for the choice maker. 
The set of options, or alternatives, a, 6, c, . . . is called the 
choice set C — {a, b,c, . . .}. The complete set of all pos- 
sible options available to any individual is called the uni- 
versal choice set and a subset of the universal choice set 
considered by a particular individual is called the reduced 
choice set. If C is the universal choice set then all pos- 
sible reduced choice sets correspond to events and form 
what is called a-algebra over C in the probability calcu- 
lus. Hereon we assume that results of a choice performed 
by any decision maker are unique for given specification 
of the choice set. This means that results of a choice 
are those elements of cr-algebra which correspond to el- 
ementary events of the probability theory. Finally, we 
consider only discrete choice models, i.e. those for which 
their choice sets arc finite and all options can be explicitly 
enumerated. 

Each alternative of the choice set can be characterized 
by a set of attributes. Generally, attributes are random 



3 



variables in the sense used in probability theory, i.e. they 
are funcions defined on reduced choice sets with values 
in other sets. We consider both directly observable at- 
tributes and functions of attributes. 

In order to construct decision rules, the decision maker 
has to be able to compare alternatives. In neocalassical 
economy Q two alternatives a and b are comparable us- 
ing the preference-indifference operator £ which orders 
C linearly, i.e. 

1. The £ is reflexive: 

V ae c a£a, (1) 

2. The £ is transitive: 

Va.f^ec (a £ b A b £ c) =► (a £ c), (2) 



3. Any two alternatives are comparable: 

Va.bec a£b y b£a. (3) 

Since the choice set C is finite, the existence of the most 
preferred alternative a* is guaranteed 

3 a - V aec a* £ a. (4) 

Because of the properties llll.'il) there exists finite random 
variable U : C — > R 1 , refered to as the utility function, 
such that 

Va,bec a £ b & U(a) > 17(6). (5) 

From eq. Q it follows that the most prefered alternative 
a* has the largest utility U 

Vaec U(a*) > U(a) (6) 

or the smallest disutility U = —U. Hence, important 
property of U 

Property 1 Disutility function is bounded from bottom. 

B. Specification of the system 

Our system S consists of a complex communication 
network Q, represented by a set of directed graphs, and 
set of decision makers. Nodes of graphs in Q represent air- 
ports and its directed edges correspond to air connections 
between them. Each ordered pair (OD) of the origin O 
and destination D ports and set of all directed paths, 
called routes, connecting O and D, excluding loops, con- 
stitute the fc-th market and is represented by the graph 
Qk (k = 1, . . . , M). The whole network Q is the sum of 
market graphs 

Q = QiU...UQm- (7) 



Using graph-theoretic terminology, markets are repre- 
sented by sets of Hamilton paths, i.e. no one site on 
the route is visited more than once. So defined markets 
can overlap. The system is further specified by defining 
all features relevant for its state as sets of potential pas- 
sengers (decision makers) for each market, resources of 
aircrafts, flight schedules, flight fares, etc. The state of 
the system is defined using results of choices performed 
by all decision makers, which means that the state-space 
consists of all possible distributions of all decision makers 
for all markets of the network. 

We assume quasistatic approximation, meaning that 
any time variation of the system as a whole is slow com- 
pared to the time of individual decision. Moreover, spec- 
ifying the state of the system, we accumulate decisions 
taken by individuals in time interval at least an order of 
magnitude shorter than the system time scale, but not 
necessarily short compared to individual decision time. 
In airline practice, one day is normally a good time bin 
for accumulation of decisions. 

In present model we do not embed the network in any 
metric space, such that only the topology, or connectivity 
scheme, is relevant and not the geometry. Accounting for 
geometry may appear necessary along with sofistication 
of our model for utility. 

For each market k we define the total market disutility 
as a sum of disutilities of all passengers from this market 
Ui k where ik stands for index of the i-th passenger in the 
fc-th market — 1, . . . , Ik) 

Uk=Y,Ui>- (8) 

ifc=l 

Overall network disutility U is a sum of market disutilities 
over all M markets 

M 

U = Y,U k - (9) 

fe=i 

States of the system are defined in discretized time, 
which means that for definitions of functions of states we 
use all individual decisions taken in given time interval, 
e.g. one day. 

C. Additivity and extensivity of the utility 

In order to find the probability of the state of the net- 
work using utility function and incorporating classical 
Boltzmann-Gibbs arguments, the utility has to be an 
additive stochastic variable. Disutility U of the system 
is called additive with respect to any two subsystems A 
and B if their disutilities add up to U, i.e. U — Ua + Ub 
(cf. e.g. ref. 8] for more detailed discussion). There- 
fore, additivity requires the notion of the subsystem to 
be clarified first. In our case, the subsystems cannot be 
defined by dividing the connectivity net only, because 
the smallest meaningful entities in the theory are the 
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markets, represented by graphs Qk (k = 1, . . . , M). We 
decompose the system S into the sum of markets 

S = Si®...®S M , (10) 

where © means the set sum of corresponding graphs 
and individuals being ascribed to corresponding markets. 
Fig. ^ illustrates decomposition of an example system 
into markets. 




. , (1,2) (1,3) CM) 

network 




(2,3) (2,4) ( 3 > 4 ) 



^ reversed (0,D) graphs 

FIG. 1: Decomposition of the network into markets 

Each market (O, D) can be further split into routes 
[O, Ci, . . . , C'l, D] in the same way. For our example 
from Fig. ^ decomposition of markets into routes is the 
following: 

(1.2) = [1,2] 

(1.3) = [1,2,3] ©[1,2,4,3] 

(1.4) = [1,2, 4] ©[1,2, 3, 4] 

(2.3) = [2, 3] ©[2, 4, 3] 

(2.4) = [2,4] © [2,3.4] 

(3,4) = [3,4] © [3,2,4] (11) 

and reversed (O, D) — ► (D, O) pairs. 

For subsystems defined in this way, from the definition 
of network disutility (|8I9[I . it follows by construction 

Property 2 Network disutility is an additive random 
function of state. 

In addition, also the extensivity of disutility, in the 
sense used in ref. |8|], can be easily justified. If n is 
the number of decision makers in the system, then ther- 
modynamic limit of disutility exists 

lira — < oo (12) 

n — >oo Ji 

for U proportional to n a (a < 1). By definition, this is 
our case with a = 1, hence 

Property 3 Network disutility is an extensive random 
function of state. 



From properties 1, 2 and 3 it follows, that network 
disutility function has all the same properties relevant 
for the construction of probability measure as the total 
energy function of finite physical system with no long- 
range interactions. Therefore it can play the same role in 
construction of statistical ensemble of systems as hamil- 
tonians do, viz. it can be used for definition of the ex- 
ponential probability measure on the state-space. For 
doing this we note that because of assumption of dis- 
aggregate choice model and additivity of disutilities, for 
any two subgroups of passengers, either belonging to dif- 
ferent markets or from the same market, disutilities of 
subgroups are independent of each other. Hence, we find 

Property 4 For any independent subgroups of decision 
makers, disutilities of subgroups are mutually indepen- 
dent random variables. 

Validity of Boltzmann and Gibbs arguments, as applied 
to U, follows from properties 1-4. Hence, the probabil- 
ity of the state of the network for given value u of the 
network disutility U is equal to 

V{U = u) ~ exp(-^it), /3>0. (13) 

An important consequence of this is that one can consider 
statistical ensemble of networks and statistical thermo- 
dynamics for them with the probability density (|13fl , and 
there is no need for non-extensive statistics with power- 
law stochastic measure. At first sight, this seems counter 
intuitive because of apparent long-range correlations in 
the communication network which usually entails usage 
of non-extensive statistics (cf. e.g. ref. 9]). This is 
understandable when realizing that only sites, and not 
markets, are really correlated in this model. Markets, 
representing elementary objects in the system, are un- 
corrected since decision makers decide on routes only, 
as mentioned above. The choice of the market is not a 
subject of the game considered here. Correlations be- 
tween markets appear when additional constraints are 
taken into account, as e.g. finite transmittance of nodes. 

In the following we restrict our model to equilib- 
rium thermodynamics, or require existence of stationary 
states, which on turn requires assumption on existence 
of the first moment of disutility function 

(U) < oo. (14) 

Otherwise, the equilibrium temperature T eq = l/(3 eq 
would not exist. We are not particularly concerned here 
neither in studying processes leading to equilibrium nor 
in all sufficient conditions for existence of equilibria in 
general. Therefore we do not start our discussion of en- 
sembles from defining transition probabilities and deriv- 
ing probability densities from detailed balance. 

In our approach we do not incorporate random inten- 
sive variables. It has been recently shown ^(| that in 
the most general case of all variables of the system be- 
ing stochastic, including intensive ones, one arrives to 
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the class of ensembles with even more general proabil- 
ity densities than Levy and Tsallis power-law functions. 
Our assumption purposefully limits generality of statis- 
tical ensembles considered here, since our original intent 
was to consider intensive variables of the system as sim- 
ple steering parameters. 



III. EQUILIBRIUM STATISTICAL ENSEMBLES 

Assuming existence of stationary regime for time evo- 
lution of complex network and using utility functions 
discussed above, it is straightforward to find probabil- 
ity measure on the state-space of the network. As known 
from statistical thermodynamics, the specific form of the 
probability density is determined by the nature of con- 
straints on the system, or the set of variables of merit 
which are allowed to be stochastic. 



A. The microcanonical ensemble 

In the simplest case of disutilities being the same for 
all states of the system or, for non-discrete case, local- 
ized in very narrow interval [Uo, Uq + SU], (5U /Uq <C 1), 
one normally assumes the probabilities for all states be- 
ing the same and calls it the principle of a priori equal 
probabilities. Ensemble of systems so defined is called the 
microcanonical ensemble and corresponds to the extremly 
tight, and rather unrealistic in our case, constraints im- 
posed on the system, and on the set of decision makers 
in particular. This means that the overall disutility is 
insensitive to the distribution of decision makers which 
could be interpreted as either lack of decision makers' 
sensitivity to the conditions of the network or as an ex- 
tremly poor offer of the carrier, giving no choice to its 
customers. This is rather trivial case and we do not dis- 
cuss it furtheron. 



B. The canonical ensemble 

As the first non-trivial case consider the canonical en- 
semble of systems where the number N of decision makers 
is fixed and the structure of connectivity, or the topology 
of the connections network, does not vary in time. One 
defines the partition function or the statistical sum of the 
system 

M 

Z{J3) = l[Z k ((3) 
fc=i 

M J k 

= iie* 



pun 



(15) 



The Zk(/3) stands for the partition function of the k- 
th market and index jk runs over all Jk distributions 



of passengers over routes belonging to the fc-th market. 
Note that each market itself is treated here as a canon- 
ical ensemble, which means that also numbers of clients 
Nk (k — 1,...,M) belonging to each market is fixed. 
This means that potential passengers from given origin 
node are at least decided as for their destination node, 
although they may not a priori know which route to fly 
there. 

Due to additivity of disutility for each market (JSJ the 
Zk can be rewritten in terms of explicit sums over Rk 
routes for the k-th market 

Z k (j3) = Z(J3,N k ) 

N k R k 

1=1 J! = l 

where ji, . . . , jN k are passenger's indices running over Rk 
routes and I runs over passengers. Assuming further the 
non- subjectivity of the utility, i.e. that utility for given 
route does not depend on the decision maker, Ui.j, = Uj 
(7 = 1,..., Nk), the partition function for the k-th market 
can be written as 



zm 



3=1 

zl(p) N * 



(17) 



where Z\ ((3) stands for the partition function for the k-th 
market with one passenger on it. 

The state of the system in the canonical ensem- 
ble is defined by specifying disutility for the complete 
set of distributions of passengers over all markets, i.e. 



{{^fi=i} 



A I 



fe=l 



Probabilities of states are 



Z(f3) 



(18) 



Noteworthly, the choice probability given by formula 118|) 
is the same as obtained in the framework of the multi- 
nomial logit choice models, usually using strong assump- 
tions on the type of distributions, extreme value or Gum- 
bel, for the stochastic components of the utility (cf. e.g. 
ref. [13). 

Mean disutility can be calculated using 1|18|) , or by dif- 
fentiation over the Lagrange multiplier (3: 

M J k 
fe=i ifc=i 

= ~^Z(fi) (19) 

and if its value (U) meas is also known from measure- 
ments, their equality determines the equilibrium temper- 
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ature j3 eq = 1/T eq . The entropy and disutility fluctua- 
tions are given by formulae 

M J k 

fe=ij h =i 

= -^ lnZ(/3) 

= (3(U) +lnZ(/3) (20) 



and 



Vor(£f) = g^lnZOT 

" -| <t7> - (21) 



C. The grand canonical ensemble 

In the grand canonical case, both disutility U and the 
number of decision makers N represent random variables, 
but the connectivity of the network remains flxed. The 
grand partition function or grand canonical sum is given 
by 



Nx,...,N M 



(22) 



where Z Nk ((3) = \~[ k=1 Z k {(3) and jl = (fit, . . . , fi M ) is the 
vector of chemical potentials for markets, corresponding 
to the vector of passengers numbers N = (Ni, . . . , Nm)- 
Here Z k is the canonical partition function for the k- 
th market (|17[) but with N k being stochastic variable. 
Following intuitive meaning of the chemical potential 
for physico-chemical systems, a priori one could think 
about different chemical potentials for different markets, 
as individual passengers from different markets may con- 
tribute unequally to the overall disutility of the system. 
Here we assume, however, the hypothesis of chemical uni- 
formity of the system which states that chemical poten- 
tials of all markets are equal: fi := fix = . . . = [Im- This 
hypothesis does not exclude multiphase systems from our 
considerations, so that our communication network cor- 
responds to the case of a one-component, though not nec- 
essarily one-phase, chemical system. In further sofistica- 
tion of this formalism one should categorize passengers 
regarding professions, wealth, aims of travels etc, and 
consider different chemical potentials for different cate- 
gories. Those characteristics are in many cases correlated 
with markets but, obviously, the scopes of those features 
are not in one-to-one relations to the markets. Formally, 
fi is the Lagrange multiplier corresponding to the con- 
straint 



(23) 



Following the same reasoning as in the canonical case, 
the grand partition function can be rewritten 

M 

= E \{^ Nk Z{P,N k ) 

N U -,N M fe=l 
M 



n 

k=l 



Z(J3, N k + l)e Nk - e-fr* 



k=l 



The state of the system in the grand canoni- 
cal ensemble is defined by specifying disutility for 
the complete set of distributions of passengers over 
all markets and for all numbers of passengers, i.e. 

{{{UjM)}£$$ =l } Nk } . Probabilities of states 



are 



Pj* wo (ft**) 



e /3(fJ.N k -U jkiNk) ) 



(25) 



This reminds again the choice probability from the multi- 
nomial logit choice model, where (3/jiN k plays the role of 
the alternative specific constant (cf. e.g. refs [fj El)- 
Here, however, this constant is specific for the whole class 
of routes for the A:-th market with N k individuals on it 
and for the particular state jk(N k ). 

First moments of U and N and their correlations can 
be found by differentiation of S over Lagrange multipliers 
(3 and fi as follows: 

M Jk(N k ) 

(p) = EE E v ]k(Nk) {p^)u 3k(Nk) 

fe=l N k Jk (N k ) = l 

= -(| in5 ^L 

M Jk(N k ) 

w = EE E PwMfiN* 

k=l N k j k (N k )=l 



Var (U) = f^lnS(/3, M ) 



(U)) 



Var(N) 



\dp 2 
= -(- 



lnS(/3, M ) 



fc=i 



(26) 
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and the entropy is given by 
S{(3,p) = - M (^1h3(Am) //3 



\d(3 (3 / /3/i 

= -^(N}+/3(U) + bxE(P,fx). (27) 

Similarily to the canonical statistics, provided the first 
moments (U) meas and (N) meas are known from measure- 
ments, the temperature 1//3 and the chemical potential 
/j, can be determined from eqns (|26(l . 



D. An outreach: possible extensions towards 
stochastic network topology — the super-canonical 
ensemble 

For the network ensembles considered so far the net- 
work topology was fixed. In real communication net- 
works it is not so and the connection graph has to be 
considered as a stochastic object. This can be realized 
by representing each pair of sites as a binary random 
variable and making the number of vertices random, and 
defining the super- canonical ensemble with new, multidi- 
mensional extensive random variable L = (L\, . . . , Lm), 
where Lk represents the total path length for the fc-th 
market (OD) 



where L = L\ + . . . + Lm will be called the total volume 
of the network. 

The state of the system in the super-canonical ensem- 
ble is defined by specifying disutility for complete set of 
distributions of passengers over all markets and for all 
possible numbers of passengers, and for all possible ran- 

dom routes, i.e. {{{^(^.^J^KjLi}^}^. 
Probabilities of states are 

e 0(fj.N k -pL k -U jk ( Nk ,L k ) ) 

V MNk , Lh) (f3,n,p) = . (32) 

The analogy with multinomial logit choice model is again 
close, with the logit alternative specific constant being 
equal to /3(/xiVfe — pLk). 

For completness, we list the formulae for first mo- 
ments and correlations of extensive variables, including 
the cummulant ((UNL)) of U, N and L accounting for 
true triple correlations between them, irreducible to pair- 
wise corrlations: 



(U) 



f—lnY(/3,fx,p) 



\d/3 
-(- 



( N ) = ^f|"lnr(/3, M ,p) 



0P-,0P 
0,P 



L(OD) 



c, 



E 



L[0,C ai ,...,C as ,D], (28) 



where cr(l), . . . , o~(S) are sequencies of permutations of 
sites corresponding to existing routes. The Lk is an anal- 
ogy of the volume and we call it the volume of the k-th 
market, represented by the graph Qk- For example, for 
the network in Fig. ^ the volumes of markets are 



Var(N) = ^-(®lnY([3,n,p) 



/3 2 \dn 2 



0,P 





= £[l,2] = 


1 








£(l,3) 


= £[1,2,3] 


+ £[1,2,4,3] 


= 2 + 


3 = 5 


cov(U,N) = 


£(1,4) 


= £[1,2,4] 


+ £[1,2,3,4] 


= 2 + 


3 = 5 




£(2,3) 


= £[2,3] + 


£[2,4,3] = 1 


+ 2 = 


3 


cov(U,L) = 


£(2,4) 


= £[2,4] + 


£[2,3,4] = 1 


+ 2 = 


3 


cov(N,L) = 


£(3,4) 


= £[3,4] + 


£[3,2,4] = 1 


+ 2 = 


3 (29) 



Var(N) = ^(JjjL.hxY<J3,»,p) 
l(d_ ( d_ 

^2 \ Qp V 



/3,M 



0,P 



lnY((3,»,p)) Mhi 



0,n 



such that the volume of the whole market, including re- 
versed (O, D) — > (D, O) pairs, amounts to 40. 
Then we define the super partition function 



and 



E 



M 

E 

N u ...,Nm k=l 



n 



-0{pL h -iiN k ) <Z k {j3) 



(30) 



((UNL)) 

= (UNL) - (UN)(L) - (UL)(N) 
-(NL)(U) + (U)(N)(L) 

KKI ( | lnF(/3 ' M ' P)) ^),^ ^ 



0,H 



where p stands for the intensive variable, coupled to ex- 
tensive variables Lk and being direct analog to the pres- 
sure, and Zk((3) is the partition function given by eq. 
I[17jl. Formally, p is the Lagrange multiplier correspond- 
ing to the constraint 



and for the entropy 



<£> = E^> 



(31) 



fc=l 



S(/3,/j,,p) 

= (3(U) - 0n(N) + pp(L) + lnY(/3, M ,p). (34) 

We note that randomization of the network topology 
includes a couple of specific issues and its complete treat- 
ment needs further clarification. The set of extensive 
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variables L is related to the network itself and not so 
much to the choice model, as it docs not a priori reflect 
any preferences of decision makers. Even statistical en- 
sembles of bare graphs, being less complex objects with 
no individuals making their choice among graph's proper- 
ties, are already complicated enough to be described with 
a few extensive variables and distinguishing interesting 
cases of the microcanonical, the canonical and the grand 
canonical, depending on the contruction procedures. The 
concept of bare random networks, unrelated to decision 
theory in the sense used here, is a subject of interest 
since pioneering works of Erdos and Renyi |T^ . Various 
procedures of randomization of graphs with certain con- 
straints on their characteristics, defining Hamilton paths 
and building probability measures for their ensembles are 
proposed e.g . in refs |13| and |14| and summarized in refs 
and lla. 



E. Escort moments 



prop 

m 



There is a useful approach to study of temperature 
dependence of thermodynamic functions, incorporating 
escort distribution |l7j | , which for the probability measure 
H is defined as 



K q 
lPi 



(i=l,...,K) 



(35) 



where pi = J A . d/j,, Aj is the z-th phase space cell from 
the set of non-overlaping cells {A;}^ covering complctly 
the whole phase space, and ggR 1 . From sum over cells 
one excludes those where pi = 0. 

Using g-parametrized escort distributions, instead of 
ordinary probability density, one can efficiently probe 
those regions of the phase space where the probability 
measure is most concentrated (large positiye q) or most 
rarified (large negative q). It was shown [l8L Il9l | that 
using Renyi entropies 



1 K 
I q = —\nY 



(36) 



instead of ordinary Kolmogorov-Shannon entropies, and 
g-derivatives |2l| 



d q F(0) _ F(qfi) - F(l3) 



d q (3 



(37) 



instead of ordinary derivatives over inverse temperature, 
one effectively evolves the partition functions and all 
thermodynamic functions. For example, for the parti- 
tion functions Z and S one gets formulae 



In 



1 



[\nZ(ql3)-qlnZ(P)} 



[In Z(q(3, qn) - q In S(/3, /x)] 



(38) 



which control evolution of thermodynamic potentials and 
response functions with temperature and chemical poten- 
tial. 

Generally, for any stochastic function of temperature 
A(/3), the temperature scaling is equivalent to modifica- 
tion of the averaging prescription 



(A(qf3)) = (A(/3)) q 
and by iterating this procedure 

(A( qi q 2 ...q n (3)) = ((... (A((3)) qn ...) q2 ) qi 



(39) 



(40) 



where (..) q stands for average over the escort distribution 
(|35|l . This is understandable because any temperature 
dependence of macroscopic observables, being themselves 
mean values or functions of them, stays in the probabil- 
ity measure. Modification of the measure p\ in case of 
Boltzmann-like pi — P(U — Ui) ~ exp(~(3ui) is just 
temperature rescaling. For q — ► 1 one recovers classi- 
cal thermodynamics. In order to follow the temperature 
evolution of ordinary mean value or variance of A with- 
out measuring it, one has to calculate corresponding q- 
moment for given temperature and vary q. Obviously, 
the same applies to any moment of the random variable 
A((3), if it exists. 

The concept of escort moments is also useful for mon- 
itoring of desired regions of the state space, e.g. those 
where experimental information is known most precisely 
or where the data exhibit the most advantageous signal- 
to-noise ratio and systematic bias is small (cf. ref. [l^ 
for detailed discussion). 



IV. APPLICATIONS 



A. The model: analytical results and numerical 
simulations 



The formalism can be directly applied for study of re- 
lations between first moments of disutilities and market 
occupancies and their correlations and their dependen- 
cies on intensive variables. For any communication net- 
work, direct numerical simulations and evaluation of par- 
tition functions and probabilities on the state space can 
be performed. In many cases, however, this may be very 
demanding computationally because of factorial depen- 
dence of numbers of configurations on the graph size. For 
particular network topologies analytical shortcuts can be 
found and also some qualitative features of the networks 
can be inferred by pure reasoning. 

We consider three network topologies (cf. Fig. [2J: 

1. The maximum connectivity network, defined by a 
complete graph, where any two sites are directly 
connected by one-segment path (Fig. 

2. The hub-and-spoke network, where every site is only 
connected to the central hub (Fig. |2b), 
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3. The spider-web, representing an intermediate case 
between the previous ones, where the hub-and- 
spoke structure is complemented by connections 
between spokes forming the ring-like structure (Fig. 

my 




FIG. 2: The maximum connectivity graph (a), the hub-and- 
spoke (b) and the spider-web (c). 

For exploratory study of basic properties of the model 
we assume the passengers disutility being proportional 
to the number of segments of the j-th route, rij (rij — 
0, 1, . . .), disregarding any additional dependencies, i.e. 
U = vq + TijV, where v is a positive constant and vq 
corresponds to disutility of making no trip and we assume 
vq = without loss of genarality. This form corresponds 
to one of the basic features of the consumers disutility, 
which is his large reluctance to choose a route with large 
number of stops and changes. Such utility function allows 
us to find some analytical results and, in addition, it 
exhibits clear physical analogy to the one-dimensional 
quantum harmonic oscillator, where the energy of its n- 
th state depends linearly on n, E n = fruj(n + ~) [n = 
0, 1, . . .). The no-trip disutility v corresponds to vacuum 
energy tkj/2. The hierarchy of states of the harmonic 
oscillator corresponds to the market consisting of non- 
forking routes of increasing lengths, as illustrated in Fig. 

This analogy is more than purely academic observation 
since the harmonic oscillator is one of the best understood 
non-trivial systems in physics and thus the market of Fig. 
0is an ideal testground for any concept developed within 
our choice model for networks. Our definition of the sys- 
tem leaves some freedom for including or not the no-trip 
alternative to the set of routes available to the choice 
maker, provided one does not consider the no-trip event 
as a double visit in given site, which would contradict 
the Hamilton property of the path. The model is per- 
haps more realistic excluding the no-trip choice (rij = 0) 
which means that the initial set of decision makers is re- 
stricted only to those determined to go. Otherwise, po- 
tential passengers are indistinguishable from those who 




FIG. 3: The hierarchy of routes of lengths n = 0, 1, . . ., corre- 
sponding to the energy spectrum of one-dimensional quantum 
harmonic osillator. Crossed n = path (tadpole) corresponds 
to exclusion of the no-trip alternative from our model. 
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FIG. 4: Mean disutility (U) vs. inverse temperature j3 (upper 
left), the v = (3fi vs. (3 (upper right), the entropy S vs. j3 
(lower left) and 5* vs. (U) (lower right) for the maximum 
connectivity graph with 5 vertices. In the lower right panel, 
the values of U = 3.062 corresponding to (3 > extend to the 
left and those for /3 < - to the right from the /3 = point 
of (U) — 3.062, indicated as vertical line. Values of entropies 
are divided by 1000. Curves correspond to analytical results 
and points are from numerical simulation. 



do not enter the game at all and for which the choice set 
C is undetermined. Closed unit-edge loops, or tadpoles, 
can be easily avoided when the network topology is de- 
fined and non-random but it is known that special care 
in ascribing probabilities to the states is needed in case 
of random topologies . 

The statistical sum is trivially calculable fo both cases 



Z((3) = 
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FIG. 5: The variance of disuility Var(U) vs. inverse tem- 
perature j3 (upper left), the variance of the number of indi- 
viduals Var(N) vs. (5 (upper right), Var(U) vs. U (middle 
left), Var(N) vs. U (middle right), the covariance of U and 
N, cov(U,N) vs. f3 (lower left) and cov(U,N) vs. U (lower 
right) for the maximum connectivity graph with 5 vertices. 
The values of (U) = 3.062 corresponding to f3 — are indi- 
cated by vertical lines for figures with U on the abcissae. 



FIG. 6: The g-mean disutility (U) q vs. q (upper left) and the 
g-mean number of individuals {N) q vs. q (upper right), the 
q- variance of disutility Var(U) q (middle left), the q- variance 
of the number of individuals Var(N) q (middle right) and the 



g-covariance of U and N, 



(U, N) (bottom left) for the 



maximum connectivity graph with 5 vertices. The ratio of 
passengers in a given type of market as a function of mean 
disutility is displayed in the lower right panel. The vertical 
line at (U) = 3.062 corresponds to j3 = 0. 



The first moments of disutility in the low- and high- 
temperature (resp. cool and hot) limits are 



(U) = 



and 



(U) 



e 0v _ 

i/7? 



i 



fiv > 1 (cool) 
^«1 (hot) 



-0v 



1/7?, 



/3v > 1 (cool) 
@v < 1 (hot) 



(42) 



(43) 



for the no-trip included and the no-trip excluded model, 
respectively. Asymptotic behaviour is the same for both 



The variance of disutility, being proportional to the 
derivative of (U), is the same for both models: 



Var{U) 



v 2 e -f3v 



(l-e-0«) a 

0, /3v > 1 (cool) 
1 1(3 2 , (3v < 1 (hot) 



(44) 



Our calculations and simulations are performed in the 
framework of grand canonical ensemble and we incorpo- 
rate the notation 



(k = l,...,M) (45) 
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Then the grand partition function i|24|) is given by 



M 
k—1 



(46) 



and formulae (|26[1 can be rewritten in more specific and 
simpler form: 



(Uk) 



(47) 



where {Uk)-\ stands for mean disutility of the one- 
passenger market k, where 



l-X, 



(48) 



From cqn gSJ it follows that X k = (N k )/((N k ) + 1) < 
1 and the grand partition function l|46(l is always non- 
singular. 

Second moments and correlations are equal to 



Var(U) = Var(U k ) 
fc=i 

M 

= J2 (( N k) Var{U k ) } + (U k ) 2 Var(N k )) 

k=l 
M 

Var{N) = E Var(N k ) 



k=l 
M 



cov(U,N) 



E + W 

k=l 

M 

^2cov(U k ,N k ) 

k=l 

M 



(49) 



k=l 



where Var(Uk)\ is the variance of disutility for the one- 
passenger market k. 

The canonical case is recovered for non-random (Nk) = 
N k for all k, when Var(N k ) = 0, as expected. 

Worthwhile to note, for non-random market distility 
Var(U k ) — 0, moments of U and N are the same as for 
the system of ideal bosons (cf. ref. 0). 



B. Case study 1: the maximum connectivity 
network 



If the network graph consists of T nodes, the total num- 
ber of markets is M — T 2 — T. Using combinatorics and 
denoting by v disutility for one segment of any route, 



1750 

<U> 1500 
1250 
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FIG. 7: Mean disutility (U) vs. inverse temperature /3 (upper 
left), the v = fin vs. f3 (upper right), the entropy S vs. ft 
(lower left) and S vs. (U) (lower right) for the hub-and-spoke 
graph with 1 hub and 4 spokes. The values of (V) = 1.630 
corresponding to f3 = are indicated by vertical lines for 
figures with U on the abcissae. 



same everywhere, one finds the partition function for the 
fc-th market 



r-i 



zm = (E 



(r-2)! 

X (T-1-0! 



-0vl 



(50) 



and Z((3) is given by eqn. (|15fl . The combinatorial factor 
accounts for the degeneration of the state of given disu- 
tility and is equal to the number of routes of the same 
length. For fixed N = Y^ =1 N k 



(u) = Nv(r-i- 



r>i 



Nv(r 



and 



Var(U) 



r>i 



Nv 2 e? v 
v{NvT - (U)) 



(51) 



(52) 



where we used equilibrium temperature (3 eq = - ln(r — 

determined from eqn. (|51|l . 
In general case we failed to find analytic formulae for 
moments of and correlations between N and U. The 
chemical potential /i can be determined, provided (N) is 
known. 

It is instructive to discuss analytic solutions for the 
simplest non-trivial case of one market on the T = 3 
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FIG. 8: The variance of disuility Var(U) vs. inverse tem- 
perature P (upper left), the variance of the number of indi- 
viduals Var(N) vs. (5 (upper right), Var(U) vs. U (middle 
left), Var(N) vs. U (middle right), the covariance of U and 
TV cov(U,N) vs. /3 (lower left) and cov(U,N) vs. U (lower 
right) for the hub-and-spoke graph with 1 hub and 4 spokes. 
The values of (U) = 1.630 corresponding to f3 = are indi- 
cated by vertical lines in figures with (U) on the abscissae. 
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FIG. 9: The q-mean disutility (U) q vs. q (upper left) and the 
g-mean number of individuals {N) q vs. q (upper right), the 
q- variance of disutility Var(U) q (middle left), the q- variance 
of the number of individuals Var(N) q (middle right) and the 
q-covariance of U and N cov q (U, N) (bottom left) for the 
hub-and-spoke graph with 1 hub and 4 spokes. The ratio of 
passengers in a given type of market as a function of mean 
disutility is displayed in the lower right panel. The vertical 
line at (U) = 1.630 corresponds to j3 = 0. 



graph and assuming v = 1 in order to elliminate triv- 
ial factors. Using eqns l|47|) and (|48|) . one easily finds 
equations for equilibrium (3 and /i 



Peg = In 



2W ~ (U) 
(U) - (N) 



(53) 



and 



= 2 In 



2(N) - {U) 



In 



(U) - (N) 



l + (N) 1 + {N) 

Using eqns H53|) . 154fl an d (|49|l . we obtain 

Var(U) = -2(N) + 3{U) + (U) 2 

= (-{U) 2 /(N) + 3{U)-2(N)) 
+ ((U) 2 /(N) + (U) 2 ) 



(54) 



(55) 



where the first term corresponds to the canonical ensem- 
ble, with fixed number of individuals (N) , and the second 
term comes from randomization of the number of individ- 
uals in the grand canonical ensemble. Since in our model 
(U) is proportional to (TV), the asymptotic behaviour of 
both terms can be investigated expanding (U) around 
(N) . Asymptotically, we find that the first term vanishes 
and the second term behaves as (/V) 2 for (TV) 2§> 1. 
In addition, we calculate for this case 



and 



Var(N) = (N)((N) + 1) 



cov(U,N) = (U)((N) + 1) 



(56) 



(57) 



and we see that both Var(N) and covariance of (U) and 
(N) depend quadratically on (N) in the limit of (N) —> 



13 



oo. For completness, the entropy in this case is equal to 



5' 



In 



1 



1 + (N) 



((N) -(U)) In 



(-2(iV) + (U}) In 

-(N) + (U) 
1+(N) 



2(N) ~ (U) 
1 + (N) 



(58) 



We also investigated analytically higher T and found 
qualitatively the same bahaviour as for T = 3. 

We performed numerical simulation for the maximum 
conectivity network with r — 5 and (N) = 1000 and 
displayed our results in Figs 0] and [SJ Both f3 and /x 
were treated formally as Lagrange multipliers and thus 
we allowed for both positive and negative values of those 
parameters. The behaviour of (U), Var(U) and S with 
temperature is intuitively appealing for /3 > 0. 

The dependence S((3) exhibits the same tendency for 
systems of bosons with inversed populations alowed, i.e. 
the entropy increases with temperature for both positive 
and nagative temperatures. Moreover, the entropy S ex- 
ists in the limits of f3 — > ±oo, where S ^ — ™ const, and 

S 13 — >°° const. ' . For the canonical ensemble const. = 
and for random N, in general const. ^ 0. 

Dependence of v on (3 exhibits steeper slope for (3 < 
which is related to less steep dependence of (U)((3) or 
(3 < 0, as compared to (3 > 0. 

It is interesting to note that Var(N) is insensitive to 
neither (3 nor (U) which reflects the maximal symmetry 
of the maximum connectivity network. Whatever the 
temperature is, relative populations for all routes are the 
same and the width of distribution of N does not depend 
on temperature. 

Results of simulations for g-moments are displayed in 
Figs El One observes the same dependence of the q- 
moments on q and explicit /3-dependence of the ordinary 
moments in Figs [S] It clearly confirms interpretation of 
q as a temperature scaling parameter, as mentioned in 
chapt. III.E. Consistently, the values of ordinary mo- 
ments for (3 — are the same as q- moments for q = 0. 

Fig. 03 (bottom right) presents the ratio of passengers 
in a given market as a function of utility, where each 
point correspons to one overall mean disutility, or tem- 
perature. Clearly, there is only one type of markets be- 
cause of the maximal symmetry of this network. Vertical 
line at (U) = 3.062 divides the domain of (17) into sub- 
regions corresponding to (3 > (to the left) and (3 < 
(to the right). 



C. Case study 2: the hub-and-spoke network 

For this topology, the one-passenger disutility in the 
A;-th market is equal to 



2v, O and D spokes (case (i)) 



otherwise 



(case (ii)) 



(59) 



Denoting by M^ 2 ) the number of markets for the case 
(i) and (ii), which for the total numbers of spokes T — 1 




FIG. 10: Mean disutility (U) vs. inverse temperature /3 (up- 
per left), the v = l3/i vs. /3 (upper right), the entropy 5* vs. 
f3 (lower left) and 5* vs. (U) (lower right) for the spider-web 
graph with 1 hub and 4 spokes. The values of (U) = 2.954 
corresponding to /3 = are indicated by vertical lines for 
figures with U on the abcissae. 



are equal to T 2 — 3T + 2 and T(r — 1), respectively, and 
using l |l 7(1 . one gets for the partition function 



Mi Mi + M 2 

-2(}vN k 77 c ~/3vN k 



(60) 



k=l 



k=Ah 



and, for the canonical cncscmblc, the mean disutility is 
equal to 



(U) = (2N + N H )v 



(61) 



where N = N k and N H = Y^llu+i are total 

numbers of individuals for cases (i) and (ii) , respectively. 
For large N^s, the grand partition function is equal to 



3(M 



1 



1 



(1 _ e /3(M-2^))Mi n_ e *-«)W 2 

and the first moments of U and N are equal to 
2vM r vM 2 



(62) 



(U) 
(N) 



e -P(H-2v) _ I e -f3(n-v) _ I 

Mi M 2 

e -/3(M"2f) _ 1 + p -f3(p-v) _ I 



(63) 



Following the same procedure as for the maximum con- 
nectivity network and assuming v = 1, one finds the vari- 
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FIG. 11: The variance of disuility Var(U) vs. inverse tem- 
perature j3 (upper left), the variance of the number of indi- 
viduals Var(N) vs. (5 (upper right), Var(U) vs. U (middle 
left), Var(N) vs. U (middle right), the covariance of U and 
N, cov(U,N) vs. f3 (lower left) and cov(U,N) vs. U (lower 
right) for the spider-web graph with 1 hub and 4 spokes. The 
values of (U) = 2.954 corresponding to j3 — are indicated 
by vertical lines for figures with U on the abcissae. 



FIG. 12: The q-mean disutility (U) q vs. q (upper left) and the 
q-mean number of individuals (N) q vs. q (upper right), the 
g-variance of disutility Var(U) q (middle left), the q-variance 
of the number of individuals Var(N) q (middle right) and the 
q-covariance of U and N cov q (U, N) (bottom left) for the 
spider-web graph with 1 hub and 4 spokes. The ratio of pas- 
sengers in a given type of market as a function of mean disu- 
tility is displayed in the lower right panel. The vertical line 
at (U) = 2.954 corresponds to j3 = 0. 



ances: 



Var(U) = 



Var (N) 



(N)(-2 + -(N)) 



(64) 



the covariance 



cov (U, N) = -(N) 2 + (U)(l --(N) + ~{U)) (65) 



and the entropy 
S = (17) In 
- (iV)ln 



(-12+(N)-(U))(2(N)-(U)) 
((N)-{U))(8 + 2(N)-(U)) 

(-12 + (N)-(U))(2(N)-(U)) 2 
((N)~(U))(8 + 2(N)-{U)) 2 



!ln- 



12 In 



8 + 2{N) - (U) 
12 



12 — (N) + (U) 



(66) 



As for the previous case, we performed numerical sim- 
ulation for the hub-and-spoke network with T = 5 (one 
hub and four spokes) and (N) = 1000 and show the 
results, together with analytic curves, in Figs [7] and |S] 
Quantitatively, the (U), v and S with f3 exhibit the same 
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type of behaviour, including asymptotics, as observed for 
the maximal graph. Dependence of (U) on (3 for the hub- 
and-spoke for small f3 is steeper than for the maximum 
connectivity graph. The maximum of entropy at f3 = 
is sharper for the hub-and-spoke than for the maximum 
connectivity graph. 

Contrary to the maximum connectivity graph, the vari- 
ance of N does depend on temperature (Fig. [S]upper and 
middle right) and decreases with temperature for (3 > 0. 
The same tendency is observed for variance of U (Fig. 
IH1 upper and middle left). This behaviour reflects lower 
symmetry of the hub-and-spoke network, as compared 
to the maximum connectivity, which means that not all 
markets are identical. There are two types of markets, 
spoke-spoke and spoke- hub (cases (i) and (ii) in eqn (JEU) 
and relative occupancy of the second one increases with 
temperature, making the market probability distribution 
narrower, eventually reaching its minimum for [3 — > 
(T — » oo). The mirror- like behaviour is observed for neg- 
ative (3 but the plateau values for Var{U), Var(N) and 
cov(U, N) are different than for positive (3. 

Simulations of q-moments are displayed in Figs |5] 
Again, simulating the g-dependence one reproduce all 
features of explicit /3-dependences of moments and cor- 
relations observed in previous figures. 

Fig. El (bottom right) presents the ratio of passengers 
in a given market as a function of utility, where each 
point corresponds to one mean overall disutility. Aster- 
isks correspond to the market with disutility v, i.e. the 
spoke- hub market, and triangles to the spoke-spoke mar- 
ket with disutility 2v. The region of (3 > extends to 
the left of the vertical line at (U) = 1.630 and j3 < - to 
the right. 



in Fig. 1131 where probability distributions for binary 
and ternary random variables are shown as functions of 
j3. For the binary case (cf. Figs El upper), the widest 




FIG. 13: Various /3-dependent probability distributions lead- 
ing to single or multiple extrema for second moments of disu- 
tilities and numbers of individuals, located at different values 
of j3. Upper panels present two types of probability distribu- 
tions for binary random variable: the maximum variance at 
13 = (upper left) and the maximum variance at (3 > (up- 
per right). Lower panels show two types of ternary random 
variables: one with maxima at /3 — and at f3 > (lower left) 
and with three local maxima, two of them for f3 > (lower 
right). 



D. Case study 3: the spider-web network 

The spider-web network represents an intermediate 
case between the maximum connectivity and the hub- 
and-spoke toplogies, discussed above. But concerning 
symmetry, the spider topology is less symmetric than pre- 
vious two, in a sense that there are three different types 
of markets in this network. We found analytic formulae 
also for this case but they are very lenghty and we do 
not print them here. All the results shown in Figs llUllTTl 
and 1 121 come from analytic formulae and from direct sim- 
ulations with the same conditions as for previous cases. 
Many qualitative characteristics of the expected values 
and g-expected values of U and N, and also of the en- 
tropy S, are the same as for more symmetric topologies, 
but inspection of Figs El and El reveals also some in- 
teresting structures in /^-dependence of the (q-)variances 
and (g-)covariances, deserving more attention. Contrary 
to more symmetric topologies, one observes two minima: 
one at (3 = or q = 0, and another one for (3 > or 
q > 0, and two maxima at non-zero values. Qualita- 
tively, such behaviour can be understood by following 
temperature dependence of probabilities, as illustrated 



distribution corresponds to equal probabilities p\ = P2 
which can happen either for (3 = or (3 > 0. More com- 
plex structure of maxima for variance can easily arise 
for multinomial probability density function where prob- 
abilities intersect pairwise in more than one regions, or 
points of intersections clusterize. An example of ternary 
random variable is presented in lower Figs ll3l where dis- 
tributions for multiple maxima are illustratively drawn. 

In order to find equilibrium values /3 eq and v eq , one 
has to solve equations for (U) and (N). We observe that 
degree of those equations depends on the symmetry of 
the network and is equal to the number of topologically 
different markets for given networks: for the maximum 
connectivity all markets are the same, for hub-and-spoke 
there are two types of markets and for the spider-web 
there are three types of markets. We further observe 
that the number of local minima is less by one than the 
degree of those equations. We strongly suspect that this 
is true for higher numbers of different markets. At the 
moment we do not have any proof for the general case 
and leave this statement as a hypothesis. 

Fig. El (bottom right) presents the ratio of passen- 
gers in a given market as a function of utility, where 
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each point corresponds to one mean overall disutility, 
or one temperature. Triangles represent market where 
O is located on the rim and D in the hub, such that 
the set of possible routes for this market can be repre- 
sented by the set of disutilities {v, 2v, 2v, 3v, 3v, 4w, 4u}. 
Asterisks correspond to the market with O and D on 
the rim and non-adjacent, such that the set of routes 
is {2v, 2v, 2v, 3v, 3v, 3v, 4w, 4w}. Finally, diamonds are 
for the market of adjacent non-hub endpoints, where 
the routes are {v, 2v, 3v, 3v, 3v, 4v, 4u, Av}. The region 
of (5 > extends to the left of the vertical line at 
(U) = 2.954 and j3 < - to the right. It is interesting 
to observe that for certain values of (U), where the hub- 
rim market occupancy increases with temperature, two 
other sorts of markets tend to decrease, reaching local 
or global maxima. Such non-trivial behaviour demon- 
strates potential power of the model developed here, for 
estimation of market occupancies. 



V. CONCLUSIONS AND FINAL REMARKS 

Summarizing, we proposed a formalism, based on the 
equilibrium statistical thermodynamics, describing com- 
plex communication systems, consisting of sets of com- 
municating nodes and individuals able to make decisions 
on how they travel. In classical approach to the problem 
of choice making, based on game theory, the key concept 
is the utility function which quantifies choice probabili- 
ties. Utility function was defined axiomatically, accord- 
ing to neoclassical economic theory We defined our sys- 
tem as a sum of markets, each market represented by a 
directed graph with two communicating endpoints and 
set of Hamilton paths connecting them. Our networks 
were not embeded in metric spaces and thus the model 
deals with topologies and not with geometry. The state 
of the system is defined by ascribing each individual to 
the market and specifying its choice of the comunication 
route. Disutility of the whole system, ascribed to each 
state of it, is defined as a sum of individual disutilities, 
integrated over all markets. We have shown that such 
disutility is an additive and extensive stochastic func- 
tion and therefore the probability measure on the state 
space is of Boltzmann type, thus reproducing the multi- 
nomial logit choice model, known from decision theory, 
using quite general and fundamental reasoning. Further, 
we considered statistical ensembles: the microcanonical, 
canonical, grand canonical and super canonical. The first 
three are commonly known in statistical mechanics and 
are rather straightforward to construct and use, whereas 
the super canonical ensemble incorporates random net- 
work topology and needs much further and more detailed 
specifications. We did not randomize intensive variables 
and kept them as nonrandom variables or control param- 
eters. Generalizations in this direction and using random 
intensive fields of temperature, pressure etc. seem quite 
natural and are certainly ahead of us in further devel- 
opments. Incorporation of geometry, i.e. using metric 



characteristics of connections like distances, shapes etc., 
are quite straightforward but not equally relevant in all 
applications. In particular, for airline networks, topology 
plays principal role, as long as intercontinental connec- 
tions are not concerned. 

Following usual procedures, known from classical equi- 
librium therodynamics, we find the entropy and first mo- 
ments of extensive random variables of the model and 
their correlations. We also propose using escort distribu- 
tions for monitoring of temperature evolution of moments 
of stochastic variables. 

We applied our model to networks with specified 
topologies and we obtained analytical results for two 
topologies: the maximum connectivity and the hub-and- 
spoke, using rather simplistic disutility function, propor- 
tional to the number of route segments. For this disutil- 
ity we found physical analogy, as the hierarchy of routes 
of lengths increasing in steps of one unit can be directly 
mapped onto the one-dimensional quantum harmonic os- 
cillator and exhibits close analogy to systems of bosons. 
For the spider-web topology the formulae can be also 
found but are very complex. 

Any network can be in principle simulated numerically. 
We performed such simulations, in the framework of the 
grand canonical ensemble for all topologies, in order to 
crosscheck our analytical findings. Our simulations, how- 
ever, were done so far for rather small networks consisting 
of five nodes and average number of thousand decision 
makers. Evaluation of statistical sums requires summing 
over all possible distributions of decision makers over all 
Hamiltonian paths on the graph and this is in general 
very demanding computationally. Since the world air 
traffic integrated over one day incorporates millions of 
passengers, its realistic simulations, even with tight con- 
straints not accounted so far, will require more efficient 
software tools. 

We see considerable economic potential of our ther- 
modynamic approach to the process of decision mak- 
ing in communicating population. After complementing 
the model with constraints making it more suitable for 
real life applications in air transportation, as bunching of 
passengers in finite-size planes, realistic flight schedules, 
availability of air space and time slots, finite node trans- 
mittance, accounting for delays and other perturbations, 
and finally, using realistic disutility function, and many 
other factors, it can be used for monitoring and fore- 
casting the state of the network. Predictions for market 
occupancies and fluctuations of numbers of individuals 
depending on passenger's disutility and all types of cor- 
relations between disutilities and occupancies arc useful 
for management of resources. 

We see three principal directions of development of our 
model: 

• Network description in the framework of super- 
canonical ensemble, i.e. using random topology, 
and incorporation of geometric characteristics of 
the network by embeding the graph in metric space, 
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• Monitoring of collective phenomena in the commu- 
nication network using entropies, response func- 
tions and an approach based on the phase space 
uniformity, developed in ref. |l8[ . 

• Transforming our choice model into the full game, 
with airline operators and carriers being also active 
players in such game. This needs finding a coupling 
mechanism between utility functions of all parties, 
being e.g. typical predator-prey coupling of the 



Lotka-Volterra type, as a good starting point for 
further developments. 
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